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Gravitational radiation from a slightly distorted black hole with ringdown waveform is well under- 
stood in general relativity. It provides a probe for direct observation of black holes and determination 
of their physical parameters, masses and angular momenta (Kerr parameters). For ringdown searches 
using data of gravitational wave detectors, matched filtering technique is useful. In this paper, we 
describe studies on problems in matched filtering analysis in realistic gravitational wave searches 
using observational data. Above all, we focus on template constructions, matches or signal-to-noise 
ratios (SNRs), detection probabilities for Galactic events, and accuracies in evaluation of waveform 
parameters or black hole hairs. In template design for matched filtering, search parameter ranges 
and template separations are determined by requirements from acceptable maximum loss of SNRs, 
detection efficiencies, and computational costs. In realistic searches using observational data, how- 
ever, effects of non-stationary noises cause decreases of SNRs, and increases of errors in waveform 
parameter determinations. These problems will potentially arise in any matched filtering searches 
for any kind of waveforms. To investigate them, we have performed matched filtering analysis for 
artificial ringdown signals which are generated with Monte-Carlo technique and injected into the 
TAMA300 observational data. We employed an efficient method to construct a bank of ringdown 
filters recently proposed by Nakano et al, and use a template bank generated from a criterion such 
that losses of SNRs of any signals do not exceed 2%. We found that this criterion is fulfilled in 
ringdown searches using TAMA300 data, by examining distribution of SNRs of simulated signals. 
It is also shown that with TAMA300 sensitivity, the detection probability for Galactic ringdown 
events is about 50% for black holes of masses greater than 20Mq with SNR > 10. The accura- 
cies in waveform parameter estimations are found to be consistent with the template spacings, and 
resolutions for black hole masses and the Kerr parameters are evaluated as a few % and ~ 40%, 
respectively. They can be improved up to < 0.9% and < 24% for events of SNR > 10 by using 
fine-meshed template bank in the hierarchical search strategy. 

PACS numbers: 95.85.Sz, 04.80Nn, 07.05.Kf, 97.60.Lf 



I. INTRODUCTION 

With recent progress in the development of gravita- 
tional wave detectors as resonant bars pj and laser in- 
terferometers 0- S CU - the detectors' sensitivities and 
operational stabilities have significantly been improved. 
Several groups have performed long-term observations 
which provided data of good quality suitable for various 
kind of astrophysical gravitational wave searches. These 
are for "chirp" signals from inspiralling compact neutron 
star binaries, burst-like signals from stellar core collapses 
in supernovae, stochastic gravitational wave background, 
and continuous waves from pulsars. Thus well orga- 
nized studies on data analysis techniques are required 
for gravitational wave detection and gravitational wave 
astronomy in near future. For some classes of gravita- 
tional wave sources, data analysis techniques have al- 
ready been investigated extensively, in particular as for 
inspiral searches. The matched filtering technique is con- 
sidered as one of the best method for this, since the grav- 
itational waveform is theoretically predicted in good ac- 



curacy. Although it is rather difficult to exploit an opti- 
mal method for gravitational wave bursts with undefined 
waveforms, several methods and filters have been pro- 
posed (e.g. iH). _ 

Gravitational radiation from a perturbed black hole, so 
called ringdown, is one of the classes of known- waveform 
gravitational wave sources. The black hole perturbation 
in the late stage of a black hole formation, initiated by 
a coalescence of a binary system for example, can be de- 
scribed by quasi-normal modes (QNMs). Therefore grav- 
itational radiation is expected as of damped sinusoidal 
waveform. This is represented in time domain by 

h(t) cx exp (~^(t - *o)) sin (27r/ c (t - t ) + fa) (1) 

where / c , Q, (f>o,to are the central frequency, the quality 
factor of the ringing, the initial phase, and the arrival 
time, respectively. The QNMs were studied extensively 
by using Regge- Wheeler- Zerilli or Teukolsky formalism 
in 1970's and 1980's (e.g. HHIIOU), and several meth- 
ods have been proposed for numerical computation of 
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them (0, 0, 0, 0| , and ^3 f° r review and refer- 
ences therein). Leaver developed a numerically stable 
method for calculation of QNMs with good accuracies 
|l6l |. and his results have been verified by many authors 
with independent schemes. For the least damped mode 
of the QNM with I = m = 2 spin- 2 spheroidal harmonics 
-2S22, Echeverria have given analytic expressions for f c 
and Q as functions of the black hole mass M and the 
non-dimensional angular momentum (the Kerr parame- 
ter) a = [0,1), as [l3| 

,/c [kHz] ~ 32(-^-) X [l-a(l-af] (2) 
Q ~ 2(1 -ay (3) 

where = 3/10,7= -9/20 and a = 63/100. These ex- 
pressions give ringdown parameters (f c , Q) with accura- 
cies of ~ 5%. For black holes with masses of 10 ~ 200A/ Q , 
the ringdown frequencies fall in the observational band 
100Hz ~ kHz of ground-based laser intcrferometric detec- 
tors. Matched filtering technique is useful to find ring- 
down signals in outputs of gravitational wave detectors 
since their waveforms are modeled in a simple form as 
equation Q in terms of the parameters {f c ,Q) which 
can be easily converted into the black hole parameters 
(M,a). 

The black hole ringdowns are very interesting in grav- 
itational wave astronomy from some aspects. First, ring- 
down detection will be a direct confirmation of the ex- 
istence of astrophysical black holes. Second, it gives an 
another "mass window" for observation of astrophysical 
objects by detecting gravitational waves. For example, 
for black hole binaries with several tens of solar-masses, 
frequencies of gravitational radiation in the inspiralling 
phase are out of the observational band of contemporary 
laser interferometers, while they can be detected via ring- 
down waves. In spite of its significance, there are fewer 
publications on ringdown searches compared to those for 
inspiral or burst searches. Creighton reported the results 
of a single filter search for a black hole with a particular 
mass and spin using data of the Caltcch 40-m interferom- 
eter E3 ■ Arnaud et al. discussed a general template con- 
struction method for matched filtering searches and its 
application to the ringdown case Recently Nakano 
et al. proposed an efficient tiling method for ringdown 
filters Sill. 

The aim of this paper is to discuss problems in matched 
filtering analysis in realistic search using observational 
data. In matched filtering with a template bank, the 
template separation is taken so that losses of signal-to- 
noise-ratios (SNRs) of any signals do not exceed a pre- 
selected value (minimal match), MM. The choice of 
MM determines detection efficiencies and waveform pa- 
rameter resolutions. Some template construction meth- 
ods for a given MM have already been proposed as EH ; 
[20I l2lj . In a realistic search, however, non-stationary 
noises in detector outputs would reduce SNRs, so it is 
not clear whether the criterion > MM is still fulfilled. 



Furthermore, the template which gives the maximum of 
the matches is not necessarily the nearest to the sig- 
nal, thus the noise effects would also appear in wave- 
form parameter resolutions. To investigate these prob- 
lems, we have carried out matched filtering searches for 
artificial ringdown signals embedded in TAMA300 data 
with Monte-Carlo technique. The ringdown signals with 
randomly selected waveform parameters are injected into 
TAMA300 data by software, and filtered with the tem- 
plate bank. We employ the template spacing method re- 
cently proposed by Nakano et al. [2Qll2l[. which provides 
an efficient tiling scheme to cover the parameter space 
with a fewer templates keeping the criterion > MM. 
For gravitational wave sources, we consider black holes 
distributed in the Galaxy, since magnitudes of ringdown 
waves from distances of Galactic scales (~ 10 kpc) are ex- 
pected enough large to be detected with contemporary in- 
terferometors as TAMA300 (see FIG.HJ. From the simu- 
lation results, we evaluated SNRs and their losses, detec- 
tion efficiencies with given SNR thresholds, and the wave- 
form parameter resolutions. We also discuss possibility 
of black hole "spectroscopy" with gravitational wave ob- 
servations by using very fine-meshed template bank and 
hierarchical search technique. The code used here is the 
same which will be used in full-scale event searches using 
TAMA300 data. The information obtained in this work 
are important for astrophysical interpretation of search 
results and to set an observational limit on the ringdown 
event rate. 

The plan of the remaining part of the paper is as fol- 
lows. In Sec. II, we briefly describe the matched fil- 
tering technique and the template construction method 
employed here. The details of our simulation studies are 
given in Sec. Ill, and the results are shown in Sec. IV. 
Section V is devoted to discussion. Finally, in Sec. VI, 
we summarize and conclude our work. 



II. MATCHED FILTERING 

A. Template spacings 

The matched filtering method is to compute a cross- 
correlation (match) between detector outputs and each of 
reference waveforms as a "template" characterized with 
waveform parameter(s), = (X 1 ^ 2 ) = (f c ,Q) in 
ringdown case, to find signals possibly embedded in the 
detector outputs. The match between a signal s(t) and 
a template h(t) with the k-th waveform parameter X£ is 
expressed in Fourier domain as 

pW) = (s\h(X»)) *Jdf Sif) ^ f f \f ] (4) 

where 5(f) and h(f) express the Fourier transforms of 
functions of time s(t) and h(t), the asterisk (*) denotes 
the complex conjugation, and 5„(|/|) is the power spec- 
trum of the detector noise. In a practical search, the 



signal s in equation is replaced by the sum with de- 
tector noise n, namely o = n + s. The template bank is 
produced so that for any signals there exists at least one 
template which gives the match greater than the mini- 
mal match, as (h(X^)\h(X^ > MM for 3 A£, where 

h is normalized as (h\h) = 1. The choice of MM is 

from statistical or astrophysical requirements (a larger 
MM reduces matches in average and hence detection ef- 
ficiencies), and computational tolerance (a smaller MM 
results in an increase of total number of templates) . Thus 
the key of template construction is how cover the param- 
eter space with a fewer templates for a given MM. 

A mismatch of two waveforms, defined as 1 — (match), 
can be understood as a distance between them. Thus to 
obtain appropriate template spacings and "tile" the pa- 
rameter space, it is required to calculate the metric of 
the space [22]. As seen from equation Q}, the metric of 
the parameter space is determined by the waveform h and 
the power spectrum of detector noise S n {\f\)- In our pre- 
vious papers (2^, El > we discussed a tiling method in the 
frequency domain, assuming the detector noise is white, 
i.e., S n (\f\) = const, in equation JJJ. In this case, the 
metric and the template tiling can be formulated analyt- 
ically, hence it is easy to program the template spacing 
scheme in a search code. 

In [20l[2l| . we also discussed the validity of a use of this 
template bank in the case of a detector's colored noise 
power spectrum. As a model of realistic noise power spec- 
trum, we used a fitting curve (without lines, structures, 
fluctuations) of the noise power spectrum of TAMA300 
at Data Taking 8 in 2003. We prepare a template bank 
using the above spacing. Signals are generated so that 
their amplitudes are normalized to unity. Then, we per- 
form the matched filtering and evaluate the maximum of 
the match for each signal. It is found that, in most cases 
(99.12% in our calculation), the losses in SNR are not 
greater than 2% which is expected from the pre-assigned 
minimal match for the white noise case. This shows that 
the template spacing can be used even in the case of col- 
ored noise power spectrum. The mean of the match is 
0.993. However, in the real data analysis, we must treat 
non-stationary noise fluctuations in detector's output. In 
the following, we discuss the simulated ringdown signals 
injected into TAMA300 data. 



B. Matched Filter Code 

The matched filtering code used here has developed 
at National Astronomical Observatory of Japan (NAOJ) 
and at Osaka City University (OCU). The TAMA300 
data are split into data chunks of about 52 sec, and 
the matched filtering code are applied for each of them. 
We use the moving-average method to obtain the noise 
power spectrum of detector outputs S n (\f\) as a fre- 
quency weighting function in equation i n the calcu- 
lation of cross-correlations j2^. We define the parame- 
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FIG. 1: TAMA300 sensitivity curves for DT6 (thin) and DT8 
(thick). Ringdown spectra are also shown for black holes of 
M = 50M o (dotted) and M = 2OOM (dashed) with a fixed 
Kerr parameter a — 0.98 and a distance d — 10 kpc. 



tcr ranges for f c and Q as 1.0 x 10 2 < f c < 2.5 x 10 3 
Hz and 2.0 < Q < 33.3, and take the minimal match 
MM = 0.98. The total number of templates M amounts 
to 628. To analyze a data chunk of length 52 sec, it takes 
about 130 seconds with the 628 filters using a single CPU 
(Intel Pcntium-IV 2.5 GHz). The TAMA300 data have 
analyzed with sixteen-node PC clusters at NAOJ and 
OCU. The time required for a full-scale search using 1000 
hours data is about one week. 



In this work, we use the template bank generated so 
that the minimal match is assured when the noise power 
spectrum is white i.e. S n (\f\) — const, in the calcula- 
tion of correlations of two waveforms by equation Q . It 
must be noted that since the noise power spectrum of a 
real detector is not stationary, the metric of the parame- 
ter space and appropriate template spacings also change 
in time. Nevertheless we apply the same template bank 
to all the data chunks used here. This is supported by 
the fact that ringdown waveforms have Fourier spectra 
of rather narrow band centered at f c (see FIG.^I, which 
is especially true for large Q ringdowns, and the match 
given by equation is expected to be less sensitive to 
whole shape of the noise power spectrum S n (\f\). This 
has been confirmed in |2ll| with model calculations using 
the TAMA300 fitted noise curve as described above (sim- 
ilar discussion is also found in 0). Thus the question 
is if the template bank based on white noise power spec- 
trum can still be valid or effective enough when applied 
to real observational data that contain time variation of 
detector's noise. 
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III. MATCHED FILTERING SIMULATION 
USING TAMA300 DATA 

A. Monte-Carlo Simulation of Galactic Ringdown 
Events 

In this section, we describe our simulations in detail. 
Randomly generated ringdown signals are injected into 
TAMA data, and filtered with the template bank. For 
gravitational wave sources, we consider Galactic black 
holes, since TAMA300 has enough sensitivity to detect 
them thus they arc suitable target for our search. The 
spatial distribution of sources in the Galaxy is assumed 
as a form in the Galacto-centric cylindrical coordinates 

dN ocexp(-2R 2 /Rl)cxp{-\z\/h z )RdRdzdip (5) 

where Rq = 4.8 kpc and h z = 1 kpc j24|. The spin axes 
of black holes are assumed to be randomly directed. It 
is rather difficult to estimate the absolute amplitudes of 
ringdown waves. Flanagan and Hughes have given an 
estimate on a fractional energy of the black hole masses 
radiated as gravitational waves to be about 3% [25| , and 
we use this value for determination of the signal ampli- 
tudes. The injection point t is randomly selected in 
time of the TAMA300 observation periods. Once to is 
given, the sky position of the source in the horizontal co- 
ordinates (0, <j)) is obtained from the source position in 
the (R, z, f) coordinates, and then the antenna pattern 
function in the direction (9, <fi) is calculated. The am- 
plitude of the signal at the detector is determined from 
the distance to the source, the antenna pattern function, 
and the angular distribution of the gravitational radia- 
tion determined by the spin- weighted spheroidal harmon- 
ics _2<S'22(cosi, a), where i is the inclination angle of the 
black hole axis seen from the observer. The ringdown 
parameters (/ c , Q) of the simulated signals are uniformly 
chosen in ranges f c = 100 ~ 2500 Hz and Q — 2 ~ 33.3 
(FIG. |21 EJl ■ The range of f c is determined considering the 
observational band of TAMA300 (see FIG.UJ). The range 
of Q corresponds to the Kerr parameter a — ~ 0.998. 
The distribution of black hole masses (mass function) of 
the simulated ringdown signals results in oc M~ 2 , since 
f c is uniformly selected and M is inversely proportional 

tO f c . 

A simulated signal is superimposed on a data chunk 
of 52 sec long at a randomly selected injection point to 
with an initial phase (FIG. 0J, and the data chunk is 
filtered with each template of the bank, and also with 
the waveform of the simulated signal itself. The inputs 
for each event are the signal parameters (/ c , Q), and the 
outputs are the matches for all of the filters used. Here 
the match for each template is given by equation us- 
ing the value of the filter output p at the local maximum 
found just around the injection point to. The SNR of an 
event is defined as the greatest in the matches, and the 
signal parameters arc determined as those of the best- 
matched template. The TAMA300 data used here is ob- 
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FIG. 2: Parameter range. We generate ringdown signals with 
parameters f c = 100 ~ 2500Hz, Q = 2 ~ 33.3 uniformly 
chosen in (f c ,Q) plane. 

tained from 2001/Sep/17 to 2001/Scp/20 in TAMA Data 
Taking 6 (DT6), and from 2003/Fcb/13 to 2003/Feb/17 
in DT8. For each data chunk, 200 simulated signals are 
examined individually, and these procedures are carried 
out for all the data chunks. The noise power spectrum of 
the detector outputs 5„(|/|) required in the calculation 
of the match is evaluated for each data chunk with the 
moving average method using the preceding data chunks. 
Note that we use the fixed template bank for all the data 
chunks, as described in section fll Al and the amplitude 
of each template is normalized so as to be (h\h) = 1 in 
equation 

This simulation method has some advantages in evalu- 
ating the detection efficiency or the detection probability 
for the Galactic events compared to the case with analyt- 
ical calculations using a Galactic modeling and a typical 
detector sensitivity. Since we use the real observational 
data, we can include the time variation of detector's noise 
fluctuation and the sensitivity. Moreover, we also follow 
the antenna direction in the Galaxy at all times, thus we 
can properly take into account of the operational condi- 
tions of the detector to any Galactic events which would 
have occurred during the observation periods. 

IV. RESULTS 

A. Detection Efficiency 

By examining the SNRs of simulated ringdown events, 
we evaluate the detection efficiency e for a given SNR 
threshold. Since we assume the sources are distributed 
in the Galaxy as equation this can be interpreted 
as the detection probability for the Galactic ringdown 
events. For example, if we detect gravitational wave 
events JV b s in observation time T and assume all of them 
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FIG. 3: A part of the template bank generated under the 
condition of minimal match MM = 0.98. The templates are 
placed along fixed Q lines (Q = 2. 6,4.4, 7.7, 13.6, 24.0) by 
using the tiling method proposed in 1201 . The oval for each 
of the templates shows the tiled region in which mismatches 
between the template and any waveforms do not exceed 2% 
in the case of white noise. There is no uncovered region in the 
parameter ranges / c = 100 ~ 2500 Hz and Q = 2 ~ 33.3. The 
ratio of the total area of the tiles to that of the parameter 
space covered is 1.57. 
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FIG. 4: Schematic view of software signal injection into real 
data (in time domain): (a) strain data n(t) from the interfer- 
ometer, (b) the strain data n(t) and a simulated signal s(t) 
with waveform parameters f c = 500 [Hz],Q = 20 injected at 
to — 100 [msec], (c) matched filter output (n + s\s). 



are of Galactic origin, the Galactic ringdown event rate 
is estimated as N Q } is / (eT). The detection probabilities 
in DT6 and DT8 as functions of the ringdown frequency 
f c are displayed in FIG. |S1 for different SNR thresholds, 
SNR > 6, > 10, and > 30. All the conditions which affect 
to the detection probabilities such as the template spac- 
ings, noise fluctuations and the position of the observer 
are properly taken into account. The detection probabil- 
ities are significantly improved in DT8, in particular for 
higher f c events. This is attributed to power-recycling 



FIG. 5: Detection probability for Galactic events, for DT6 
(open symbols) and DT8 (closed symbols). Circle: SNR > 6, 
Square: SNR > 10, Triangle: SNR > 30 

implemented after DT6 to increase the input laser power 
in the cavities of the interferometer and reduce the shot 
noise which is dominant in higher frequencies. We can 
detect about e = 50% of the Galactic ringdown events of 
frequencies f c < 600 Hz (corresponding black hole mass 
> 20 M Q ) under the criterion SNR > 10. 

We note that these results depend on the evaluation of 
absolute amplitudes of the ringdown waves at the sources. 
We assume that the fractional mass loss of a black hole by 
radiating gravitational waves is 3% according to Flana- 
gan and Hughes, and all of the energies are carried by 
the waves of the least damped mode of the QNMs. 



B. Matched Filtering and Parameter 
Determination 

Next we consider matches between waveforms and tem- 
plates, and accuracies in determination of the waveform 
parameters. It is expected that these arc primarily lim- 
ited by the template spacings. We recall here the defini- 
tion of the match between a template h and a signal s 
injected in detector output n: 

P=<o\h)«Jtf lmr (6) 

where 5(f) = h(f) + §(f). For each of our simulated sig- 
nals, we can define three points in the parameter space, 
those arc (1) a priori template with the same param- 
eter values of the simulated signal, S : (fc ,Q S 'i P S = 
(o|§)), (2) the nearest to S in the template bank, N : 
(fc,Q N ;p N = (o|N», and (3) the best-matched tem- 
plate, M : (f™,Q M ;p M = (o\M)), which is the maxi- 
mum of the match in the bank and gives the SNR for the 
signal. The template we can single out from the bank 
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is M, not the nearest N, but the implicit assumption in 
matched filtering is that the two templates N, M are iden- 
tical and the parameters = f^,Q N = Q M are the 
best estimates for S. Note that the metric of the param- 
eter space which defines the distances between the two 
templates depends on the noise power spectrum S n (\f\)- 
In the template construction we employ in this work, the 
metric of the parameter space is given on the assumption 
of white noise, 5„(|/|) = constant. The template spacing 
is carried out in order that for any signals of waveform 
parameters within the search range, we can find at least 
one template such that the match is greater than the 
minimal match MM, i.e. p M / p s > MM. Thus if the 
noise power spectrum is colored, the definition of the dis- 
tance is changed and hence the maximum of the match 
p M I ' p s could be smaller than MM. Nakano et al. has 
investigated the matches between the templates and sig- 
nals with randomly selected waveform parameters, (s|N), 
by direct integration of equation JJjJ using an analytical 
fit of the TAMA300 typical noise curve (no lines, no fine 
structures), and showed that the fraction of signals which 
give p M /p s < MM is only 0.88% ([23, see also section 
ITTa) . Therefore we use the fixed template bank in the 
matched filtering analysis using TAMA300 data, for all 
the data chunks at different detector conditions. 

However, the noise power spectrum of a real detector 
S n (\f\) has many lines and fine structures, and changes 
its shape with time. Moreover, the detector noises are not 
stationary and signals would be contaminated by burst- 
like noises. Thus it is not clear whether the criterion 
p M I p > MAI is still fulfilled in a practical search using 
the fixed template bank. An alternative method is use 
of adaptive template banks, determined from the actual 
noise power spectrum and metric of the parameter space, 
as performed in neutron star binary search |26j . Our 
approach makes the data analysis codes simpler, but the 
effect of noise fluctuations must be considered carefully. 
Furthermore, when the noise power spectrum <S ra (|/|) is 
deformed and the signal s(/) in equation JHJ falls in a 
noisy frequency region, the match is degraded and the 
template M selected from the bank does not necessarily 
coincide with the nearest template N (see FIG. EJ). In 
this case, the errors in parameter determination could be 
larger than those limited by the template spacings. 



1. Finding Best-Matched Template 

To investigate these effects, we examine the ratios 
p N /p s = (o\N) I (o|S) and p M / p s = (o\M) / (o|S), those 
are measures of the utility of the template bank. If the 
noise power spectrum is white (<SVi(|/|) = const.) and 
there are no non-stationary noises, the maximum value 
of p N I ' p s is unity and its distribution has a cut-off at the 
minimal match MM = 0.98 defined in the template con- 
struction. The distribution of p N / ' p s for the TAMA300 
data-injected signals are shown in FIG. 0(a). For a ref- 
erence the distribution of (s\N) calculated by using an 



analytical fit of the TAMA300 typical noise curve [2l| is 
also displayed. The tail of the distribution of p N / p s in 
the lower side looks in agreement to that expected from 
the template design, (s|N). The fraction of events with 
p N /p s < MM is 3.3%, which is larger than 0.88% for 
(s\N) because of the noise effects. We note that about 
10% of simulated signals have the fractional match p N / p s 
greater than 1 (many of them are observed in lower Q 
signals), and the mean of p N j ' p s is 0.994 , slightly larger 
than that of (s|N), 0.993. This is related to contami- 
nations of the signals by the noise n, and the estimated 
noise power spectrum S n (\f\) used in matched filtering. 
If the true noise power spectrum at the time when the 
data n was obtained differs from the estimated spectrum 
^n(|/|)j P N I P S > 1 is possible since each template is 
normalized as (h\h) — 1 by using the estimated spec- 
trum S n (\f\). Therefore in a realistic gravitational wave 
search there are inevitable ambiguities in evaluation of 
matches (or SNRs) and their losses. 

In FIG. (b), the distribution of the ratio p M /p s is 
shown. The mean of p M / p s is 0.998, and the fractions of 
events of p M / p s < MM and p M /p s > 1 are 0.0012% and 
24%, respectively. This indicates that our template bank 
fulfills the requirement of the minimal match > MM 
even in practical search, and also it is difficult to estimate 
the "true" value of the match, (s\S) or (s|N), within accu- 
racies of a few %. The occurrence probability of N ^ M 
is about 30% (see FIG.|5J). This causes increases of wave- 
form parameter determination errors. 

2. Waveform Parameter Errors 

The accuracies in determination of the waveform 
parameters (/ c , Q) are evaluated by comparing those 
of the simulated signal § and of the best-matched 
template M. The root mean squares (RMSs) of 
the errors (A/ C // C ) RMS = ((/ C M ~ / C S )// C S ) RMS and 
(AQ/Q) RMS = ({Q M - Q S )/Q S ) RMS are calculated as 
1.3% and 22%, respectively. These are comparable to 
RMSs of (A/ C // C ) RMS = 1.22% and (AQ/Q) RMS = 
16.3% estimated from the template spacings, but some- 
what broadened. This is because for not all the events 
the template M coincides with the nearest N, as seen in 
FIG HI In practice, (A/ C // C ) RMS and (AQ/Q) RMS de- 
pend on Q, as the template distances are (see FIG. [3J ■ 
The values of (A/ C // C ) RMS and (AQ/Q) RMS are listed 
in TABLE HI 

V. DISCUSSION 

A. Black Hole Spectroscopy 

The distances between the ringdown filters determined 
from the requirement to the minimal match MM limit 
the waveform parameter resolutions. The resolutions for 
the ringdown parameters in matched filtering analysis 
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FIG. 6: An example of contour map of SNRs in the param- 
eter space. The star is the simulated signal §, which falls 
at (f c — 1197.6Hz, Q = 9.6), and the solid points represent 
the templates around S. The circle indicates the nearest tem- 
plate N, and the square shows the best-matched template M 
which gives the maximum of SNR. In most of events the two 
templates N and M coincide, but differ in this sample. This 
is due to an additional "curvature" of the parameter space 
introduced by the noise power spectrum S„(\f\), and con- 
tamination of the signal by non-stationary noises. 




FIG. 7: Distributions of 

„M I 8 _ 



ratios of the SNRs, p N /p s = 
(o|N) /<o|§) and p M /p* = (o|M) / (o|S) (see text for defini- 
tions). The broken line is distribution of (s|N), where s is a 
randomly generated ringdown waveform, obtained by direct 
integration of equation @ using a typical TAMA300 noise 
curve 12111. 



with the template bank of MM = 0.98 using TAMA300 



data are (A/ c // c ) 



RMS 



1.3% and (AQ/Q) 



i:ms 



22% 

as described in the previous section. The question here is 
how accurately we can determine the source parameters, 
the mass M and the Kerr parameter a, from the results 
of matched filtering analysis. By using the Echeverria's 



FIG. 8: Distribution of the templates which give the maxi- 
mum of SNRs in the bank. 



RMS 
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expressions (J2J and JSJ, we found (AM/M) 
20% for low Q signals, and a few % for higher Q ring- 
downs. The error in the Kerr parameter determination 
is about (Aa/a) RMS = 40 ~ 60%. The calculated values 
of (AM/M) RMg and (Aa/a) RMS are listed in TABLE [I] 
Obviously these results depend on the template de- 
sign, i.e. template spacings. It is expected that the 
parameter resolution can be improved as template dis- 
tances become closer, and saturates when the mismatches 
between signals and neighboring templates approach to 
intrinsic ambiguities limited by the noise fluctuations. 
In case of Gaussian noise, the intrinsic parameter esti- 
mation errors can be evaluated by using the Fisher in- 
formation matrix |27l l28j| . For example, for ringdown 
signals of (f c = 580Hz, Q = 24), these are calculated 
as functions of SNRs, as (A/ C // C ) RMS = 0.30% and 
(AQ/Q) 



RMS 



14% at SNR = 10 and inversely pro- 
portional to SNRs (see Appendix B. in |2l|). To confirm 
this, we performed matched filtering simulations for fixed 
parameter ringdown signals at f c = 580 Hz, Q = 24.0 
(M = 50M Q ,a = 0.996), with very fine-meshed tem- 
plates uniformly distributed around the signal. We found 



(A/c//c) 



RMS 



0.21% and (AQ/Q) 



RMS 



9.3% for 

signals of SNR = 10-20 (FIG EjT" The parame- 
ter estimation errors are in agreement with the predic- 
tions of the Fisher information matrix, and small ex- 
cursions from them are considered due to non-Gaussian 
component of the TAMA300 data. With these results, 
the black hole parameter estimation errors arc evalu- 



ated as (AM/M) RMS 
SNR - 10, and (AM/M) 



0.9%, (Aa/a) 



RMS 



RMS 



0.3%, (Aa/a) 



24% for 
7% 



RMS 



for SNR ~ 50. No further improvements are expected 
even in a case of more finer mesh, since these limits are 
from the intrinsic ambiguities in matches due to the noise 
fluctuations. (Note that since the accuracies in the con- 
versions from (/ C ,Q) to (M, a) are only ~ 5%, we need 
more accurate formulae.) 
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Q M 


(A/c//c) RMfi 


(AQ/Q) HMS 


(AM/M) RMS 


(Ao/a) RMS 


All 


1.3 (1.2) 


22 (16) 






2.55 


8.1 (2.6) 


22 (16) 


22 (12) 


64 (35) 


4.41 


4.0 (1.6) 


24 (16) 


13 (6.6) 


41(35) 


7.70 


1.6 (1.0) 


21 (16) 


6.8 (3.9) 


39 (36) 


13.6 


0.77 (0.58) 


19 (16) 


3.1 (2.4) 


40 (36) 


24.0 


0.39 (0.33) 


19 (17) 


1.9 (1.6) 


41 (37) 



TABLE I: RMS errors in estimations of ringdown - black hole 
parameters (in %) by using the template bank of MM = 
0.98. The quantities are defined as AX/X = (X M - 
X s ) /X s , (X = f c ,Q,M,a). For comparison, analytic esti- 
mations from the template distances, which should be ob- 
tained in the absence of noise fluctuations, are also shown 
(parenthesized). The parameter conversions (fc,Q) *-> (M,a) 
are based on the Echeverria's expressions in equations (3J and 
0. 



When wc find event candidates of ringdown signals, we 
can determine the black hole parameters within the accu- 
racies described above by setting a fine-meshed template 
bank locally at each of the event candidates. There is no 
need to cover whole parameter space with the template 
bank of such closer distances, the hierarchical search 
strategy copes with computational costs j2^. The res- 
olutions for the black hole parameters obtained here are 
remarkably good, compared to those estimated dynam- 
ically for known black hole candidates in X-ray bina- 
ries or radio sources, typically 10 ~ 100% [3(j ■ Some 
X-ray binaries as GRS1915+105 eject supcrluminal jet, 
and rotating black holes are considered as the central 
engine |3ll |. Precise determination of the Kerr param- 
eter a of black holes by observing gravitational waves 
will giv e strong constrains on jet formation mechanisms 
[32l I33L l34j . Thus gravitational wave observations pro- 
vide not only a probe for black hole detections, but also 
a promising clue to establish black hole "spectroscopy" . 



VI. CONCLUSIONS 

Black hole ringdowns are promising source of gravi- 
tational radiation, and the matched filtering technique 
is suitable to extract ringdown signals from outputs of 
gravitational wave detectors. In this paper, we studied 
template constructions, detection efficiencies and wave- 
form parameter estimations in matched filter searches 
using TAMA300 data. In template construction, we em- 
ployed a simple and efficient tiling method proposed by 
Nakano et al. [20I l2lj . The template bank constructed 
with this method assures the minimal match criterion 
for any ringdown waveforms in the cases of a flat noise 
power spectrum S n (\f\) = const, and a colored spectrum 
fitted to TAMA300 typical noise curve H3, EH- The er- 
rors in waveform parameter estimation are determined 
by the template separations hence the choice of MM. In 
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FIG. 9: Waveform parameter resolutions in the case of a finely 
meshed template bank (closed circle). The simulation is car- 
ried out for a fixed waveform with parameters f c — 580Hz 
and Q = 24. The templates are uniformly placed around the 
signal with separations of A/ c = 0.13Hz (A/ c // c = 0.023%) 
and AQ = 0.2(AQ/Q = 0.83%) (broken line). The predic- 
tions from the Fisher information matrix, which are expected 
in cases of Gaussian noises are also shown (solid line). 



gravitational wave searches using real observational data, 
however, it is not clear whether the minimal match cri- 
terion is fulfilled since non-stationary noises contaminate 
signals and the time-varying noise power spectra change 
matches or distances between signals and the templates. 
In a such case, the nearest template and the best-matched 
template for a signal do not necessarily coincide and pa- 
rameter estimation errors increases. 

To investigate these effects, we performed matched fil- 
tering analysis for artificial ringdown signals injected into 
TAMA300 data with Monte-Carlo technique. We con- 
sider black holes with masses 10 ~ 200A/ Q distributed 
in the Galaxy as "realistic" search targets. With TAMA 
DT8 sensitivity and using our template bank, the de- 
tection probability is evaluated as ~ 50% for Galactic 
ringdown events of f c < 600Hz, which correspond to 
black holes with masses > 20M©, assuming the fractional 
black holes mass ~ 3% radiated as gravitational waves. 
We confirmed that the loss of SNR is less than 2% in 
most cases for ringdown signals embedded in TAMA300 
data, therefore we can use the template bank in gravita- 
tional wave searches using real observational data. We 
also showed that it is difficult to evaluate matches p or 
SNRs within an accuracy of a few % because of the sig- 
nal contaminations and the difficulties in estimation of 
"true" noise power spectrum at given observation time. 
The RMS errors in the waveform parameter estimations 
are examined by comparing the simulated signals and 
the templates singled out from the bank. We obtained 
(A/ c //c) RMS = 1-3% and (AQ/Q) RMS = 22%, which 
arc comparable with the expectations from the template 
spacings, Af c /f c = 1.22% and AQ/Q — 16.3%, respec- 
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tively. The informations obtained in this work as the de- 
tection probability and the parameter estimation errors 
will be used in statistical and astrophysical interpreta- 
tion of results from the TAMA300 full-scale event search, 
which is under progress. The key in the event search is 
to develop efficient vetoing techniques to distinguish real 
signals and spurious triggers due to impulsive noises |35l| . 

Our interest extends to how accurately we can de- 
termine the waveform parameters hence the black hole 
masses and the angular momenta by using an extremely 
fine-meshed template bank. For ringdown signals of 
f c = 580 Hz, Q = 24.0, the ultimate resolutions are eval- 
uated as (A/ C // C ) RMS = 0.21% and (AQ/Q) RMS = 9.3% 
for signals of SNR = 10 ~ 20, which are consistent 
with the predictions of the Fisher information matrix. 
The errors in black hole parameter estimations are cal- 
culated as (AM/M) RMS ~ 0.9%, (Aa/a) RMS - 24% for 



SNR ~ 10, and (AM/M) 



RMS 



0.3%, (Aa/a) 



RMS 



7% 

for SNR ~ 50. In a practical search, we will employ the 
hierarchical search strategy, performing the first event 
search with an "ordinary" template bank, of minimal 
match MM = 0.98 for example, and the second search 
for event candidates with a locally fine-meshed template 
bank. The remarkably good resolutions for (M, a) are 
heartening to understand gravitational physics and as- 
trophysics of black holes. 



emitted in the merger phase, followed by the ringdown 
signal in the last stage. To investigate effects of the pre- 
ceding waves, we performed a test of matched filtering for 
a such composite waveform injected in TAMA300 data. 
The aim of this is not to generate realistic gravitational 
waveform from a binary coalescence, nor to use composite 
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FIG. 10: An example of ringdown filter outputs in time series 
for a ringdown signal (solid) and for a composite waveform 
(broken) injected in a TAMA300 data chunk. 
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APPENDIX A: PRECEDING WAVEFORMS 

So far we have assumed that there are no preceding 
waves before the ringdown signal arrives. If the ring- 
down is initiated by a binary coalescence, however, we 
would receive gravitational waves of characteristic wave- 
forms successively, as the chirp signal at first, the waves 



waveforms as templates, but to examine whether we can 
extract ringdown signals alone from the composite wave- 
form embedded in detector outputs by using our ring- 
down filters. Here we considered binary coalescence of 
masses 25M© — 25M Q to form a black hole of 50M Q with 
a Kerr parameter a = 0.98 (f c = 520Hz, Q = 12). In this 
case, the cut-off frequency of the inspiral phase is below 
100Hz, therefore the contribution of the chirp signal is 
small. Although the waveforms in the merger phase are 
quite uncertain, we presume sine-Gaussian waveform, as 
commonly used in the burst search studies, of the cen- 
tral frequency ~ 300Hz with amplitude of twice of the 
ringdown signal and a phase to be connected analyti- 
cally with the following ringdown. We found that (i) the 
maximum of SNR is found at the template nearest to the 
simulated signal, (ii) the peak of the filter output in time 
series is found at which ringdown begins, and (iii) the 
value of SNR is consistent within 1% compared to that 
in the case of no preceding waveforms injected (FIG. ITU)l . 
This is just a preliminary inspection on effects of the pre- 
ceding waveforms in a ringdown search, and we need to 
know more physics of the merger phase for quantitative 
examinations. 
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